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Abstract 

Different change-point type models encountered in statistical in- 
ference for stochastic processes give rise to different limiting likelihood 
ratio processes. In this paper we consider two such likelihood ratios. 
The first one is an exponential functional of a two-sided Poisson pro- 
cess driven by some parameter, while the second one is an exponential 
functional of a two-sided Brownian motion. We establish that for suffi- 
ciently small values of the parameter, the Poisson type likelihood ratio 
can be approximated by the Brownian type one. As a consequence, 
several statistically interesting quantities (such as limiting variances 
of different estimators) related to the first likelihood ratio can also be 
approximated by those related to the second one. Finally, we discuss 
the asymptotics of the large values of the parameter and illustrate the 
results by numerical simulations. 
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1 Introduction 



Different change-point type models encountered in statistical inference for 
stochastic processes give rise to different limiting likelihood ratio processes. 
In this paper we consider two of these processes. The first one is the random 
process Z p on IR defined by 

{p U+(x) — x, if x > 0, 

(1) 
— pll_( — x) — x, if x ^ 0, 

where p > 0, and n + and Il_ are two independent Poisson processes on R + 
with intensities l/(e p — 1) and 1/(1 — e~ p ) respectively. We also consider the 
random variables 

x Z p {x ) dx 

C P = r ry ( s , and £p = argsupZ p (x) (2) 
J K Z p {x) dx xeR 



related to this process, as well as to their second moments B p = and 

m p = Ee p . 

The process Z p (up to a linear time change) arises in some non-regular, 
namely change-point type, statistical models as the limiting likelihood ratio 
process, and the variables ( p and £ p (up to a multiplicative constant) as the 
limiting distributions of the Bayesian estimators and of the maximum like- 
lihood estimator respectively. In particular, B p and M p (up to the square 
of the above multiplicative constant) are the limiting variances of these esti- 
mators, and the Bayesian estimators being asymptotically efficient, the ratio 
E p = B p /M p is the asymptotic efficiency of the maximum likelihood estima- 
tor in these models. 

The main such model is the below detailed model of i.i.d. observations in 
the situation when their density has a jump (is discontinuous). Probably the 
first general result about this model goes back to Chernoff and Rubin pQ. 
Later, it was exhaustively studied by Ibragimov and Khasminskii in [TOl 
Chapter 5] (see also their previous works [7j and [8]). 



Model 1. Consider the problem of estimation of the location parameter 9 
based on the observation X n = (X±, . . . ,X n ) of the i.i.d. sample from the 
density f{x — 9), where the known function / is smooth enough everywhere 
except at 0, and in we have 

^ lim f(x) = a^b = \im f(x) ^ 0. 
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Denote P$ the distribution (corresponding to the parameter 9) of the 
observation X n . As n — > oo, the normalized likelihood ratio process of this 
model defined by 

dP^.u " f (X - 9 - ±) 

converges weakly in the space Po( — oo, +00) (the Skorohod space of functions 
on IR without discontinuities of the second kind and vanishing at infinity) to 
the process Z a ^ on IR defined by 



\nZ a fi(u) 



m(f)n 6 (u) - (a-b)u, ifu^O, 
- ln(f ) n o (-u) - (a - 6) u, if u ^ 0, 



where Il b and Il a are two independent Poisson processes on IR + with intensi- 
ties b and a respectively. The limiting distributions of the Bayesian estimators 
and of the maximum likelihood estimator are given by 

L u Z a b(u) du 
^ b = r ry '/ x , and ia,b = argsupZ aj6 (M) 
J R Z afi [u) du uGR 

respectively. The convergence of moments also holds, and the Bayesian es- 
timators are asymptotically efficient. So, E^ 6 and E£^ b are the limiting 
variances of these estimators, and E£^ b /E£^ 6 is the asymptotic efficiency of 
the maximum likelihood estimator. 

Now let us note, that up to a linear time change, the process Z a ^ is 
nothing but the process Z p with p = |m(|)|. Indeed, by putting u = ^ we 
get 



h\Z afi (u) 



- ln (f)n a (-^)-a;, if ^5 <0, 
In Z p (x) = In Zpi^a — b)uj. 



So, we have 



and hence 



Ca,b = ~^~~T aIld ta,b = ~^~T > 

a — a — o 



Some other models where the process Z p arises occur in the statistical 
inference for inhomogeneous Poisson processes, in the situation when their 
intensity function has a jump (is discontinuous). In Kutoyants [HJ Chapter 5] 
(see also his previous work |12] j one can find several examples, one of which 
is detailed below. 



Model 2. Consider the problem of estimation of the location parameter 
9 G ]a,/3[, < a < (3 < t, based on the observation X T on [0, T] of the 
Poisson process with r-periodic strictly positive intensity function S(t + 9), 
where the known function S is smooth enough everywhere except at points 
t* + rk, k G Z, with some t* G [0, r], in which we have 

4 limS(t) = 5_ 4 S + = limSYt) ^ 0. 

Denote Pj the distribution (corresponding to the parameter 9) of the 
observation X T . As T —>■ oo, the normalized likelihood ratio process of this 
model defined by 

Zt{u) = ~wf {x ] = exp u ln ^r dx{ty I ^ (t) ~ Se{t) ^ dt 

converges weakly in the space T> (— oo, +oo) to the process Z Tt s_ ) s + on R 
defined by 



InZ 



r,S_,S 4 



ln(^)n 5 _(a)-(5 + -5_)a, ifu>0, 



where Us_ and II s + are two independent Poisson processes on R + with in- 
tensities 5_ and S + respectively. The limiting distributions of the Bayesian 
estimators and of the maximum likelihood estimator are given by 

(r,s-,s+ = —r~^ / \ j and £r,s-,s+ = argsupZ TiS _ iS (u) 

J R Z Tt s„,s + {u) du ue u 

respectively. The convergence of moments also holds, and the Bayesian es- 
timators are asymptotically efficient. So, ECt,S-,s+ and ^£t,S-,s+ are ^ ne 
limiting variances of these estimators, and E£^ s s /E£^ s s is the asymp- 
totic efficiency of the maximum likelihood estimator. 
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Now let us note, that up to a linear time change, the process Z T: s_,s+ is 

u 



nothing but the process Z p with p = | In ( ^ ) |. Indeed, by putting 
we get 

Z r ,s-,s+(u) = Z p (x) = Z p ( S+ ^_ S ~ u\ 



So, we have 



and hence 



(r,S-,S+ - g _ g all< ^ (r,S-,S+ $ - S 



EC r ,5_,5 + - (fif+ _ ^_ )2 , EC TiS _, s+ - {s+ _ s _ )2 and ^— - _ ^ . 



The second limiting likelihood ratio process considered in this paper is 
the random process 

Z (ar) = exp jw(a;) - ^ , x e R, (3) 

where is a standard two-sided Brownian motion. In this case, the limit- 
ing distributions of the Bayesian estimators and of the maximum likelihood 
estimator (up to a multiplicative constant) are given by 

frr, X Zq^X^ djX 

Co = r ry , s , and £ = argsupZ (x) (4) 
J K Z (x) dx xeK 

respectively, and the limiting variances of these estimators (up to the square 
of the above multiplicative constant) are B = E£q and M = E£q. 

The models where the process Zq arises occur in various fields of statis- 
tical inference for stochastic processes. A well-known example is the below 
detailed model of a discontinuous signal in a white Gaussian noise exhaus- 
tively studied by Ibragimov and Khasminskii in [TUJ Chapter 7.2] (see also 
their previous work |9j), but one can also cite change-point type models of 
dynamical systems with small noise (see Kutoyants p2] and p31 Chapter 5]), 
those of ergodic diffusion processes (see Kutoyants [151 Chapter 3]) , a change- 
point type model of delay equations (see Kiichler and Kutoyants [H]), an 
i.i.d. change-point type model (see Deshayes and Picard [3]), a model of a 
discontinuous periodic signal in a time inhomogeneous diffusion (see Hopfner 
and Kutoyants [6]), and so on. 
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Model 3. Consider the problem of estimation of the location parameter 
9 E ]a,j3[, < a < P < 1, based on the observation X £ on [0,1] of the 
random process satisfying the stochastic differential equation 

dX £ (t) = -S(t- 9) dt + dW(t), 

where W is a standard Brownian motion, and S is a known function having 
a bounded derivative on ]— 1,0[U]0,1[ and satisfying 

\imS(t) - \im S(t) = r ^ 0. 

Denote P| the distribution (corresponding to the parameter 9) of the 
observation X s . As e — > 0, the normalized likelihood ratio process of this 
model defined by 

dP £ 

Ze{u) = 

= expji J [S(t -9- e 2 u) - S(t - 9)] dW{t) 

[ s ( t - 9 - £2 ^~ s ( t - e rf dt } 

converges weakly in the space Cq(— oo, +oo) (the space of continuous func- 
tions vanishing at infinity equipped with the supremum norm) to the process 
Zq(t 2 u), u G R. The limiting distributions of the Bayesian estimators and 
of the maximum likelihood estimator are r~ 2 ( and r~ 2 £ respectively. The 
convergence of moments also holds, and the Bayesian estimators are asymp- 
totically efficient. So, r~ 4 Bo and r _4 Mo are the limiting variances of these 
estimators, and Eq is the asymptotic efficiency of the maximum likelihood 
estimator. 



Let us also note that Terent'yev in [20] determined explicitly the distri- 
bution of £ and calculated the constant M = 26. These results were taken 
up by Ibragimov and Khasminskii in [TD1 Chapter 7.3], where by means 
of numerical simulation they equally showed that Bq = 19.5 ± 0.5, and so 
Eq = 0.73 ± 0.03. Later in [5], Golubev expressed B in terms of the sec- 
ond derivative (with respect to a parameter) of an improper integral of a 
composite function of modified Hankel and Bessel functions. Finally in [18], 
Rubin and Song obtained the exact values B = 16 C(3) an d E = 8£(3)/13, 
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where ( is Riemann's zeta function defined by 

oo 1 
n=l 

The random variables C p and £ p and the quantities B p , M p and £7 p , p > 0, 
are much less studied. One can cite Pflug PS] for some results about the 
distribution of the random variables 

argsupZ p (a;) and argsupZ p (a;) 

related to £ p . 

In this paper we establish that the limiting likelihood ratio processes 
Z p and Z are related. More precisely, we show that as p — > 0, the pro- 
cess Z p (y/p), y G R, converges weakly in the space V (—oc, +00) to the 
process Z Q . So, the random variables p( p and p£ p converge weakly to the 
random variables Co an d £0 respectively. We show equally that the conver- 
gence of moments of these random variables holds, that is, p 2 B p — > 16 £(3), 
p 2 M p -> 26 and E p -> 8 C(3)/13. 

These are the main results of the present paper, and they are presented 
in Section [21 where we also briefly discuss the second possible asymptotics 
p — > +00. The necessary lemmas are proved in Section [3j Finally, some 
numerical simulations of the quantities B p , M p and E p for p G ]0,oo[ are 
presented in Section HI 



2 Main results 

Consider the process X p (y) = Z p (y/p), y G R, where p > and Z p is defined 
by ©. Note that 

f y- 7 x , = P Cp and argsup X p (y ) = p £ p , 

where the random variables ( p and £ p are defined by (j2j). Remind also the 
process Z on IR defined by (EJ) and the random variables Co and £0 defined 
by dU). Recall finally the quantities 5 P = E<C p , M p = Eg, E p = B p /M p , 
B = EC0 2 = 16C(3), M = E£2 = 26 and £ = B /M = 8C(3)/13. Now we 
can state the main result of the present paper. 

Theorem 1 The process X p converges weakly in the space T>o(—oo, +00) to 
the process Zq as p — > 0. In particular, the random variables pC, p and p£ p 
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converge weakly to the random variables (q and £o respectively. Moreover, for 
any k > we have 

p k E( k p ^E( k and ^E#->E& 
and in particular p 2 B p — ► 16^(3), p 2 M p — > 26 and — > 8£(3)/13. 

The results concerning the random variable ( p are direct consequence 
of Ibragimov and Khasminskii [TOj Theorem 1.10.2] and the following three 
lemmas. 

Lemma 2 The finite- dimensional distributions of the process X p converge 
to those of Zq as p — > 0. 

Lemma 3 For all p > and all j/^^eRwe have 

^Xl^-Xl'^y^ ^\\y x -y 2 \. 

Lemma 4 For any c £ ] 0, 1/8 [ we have 

EX p 1 / 2 (y)^exp(-c| 2 /|) 
for all sufficiently small p and ally 6 R. 

Note that these lemmas are not sufficient to establish the weak conver- 
gence of the process X p in the space V Q (— oo, +oo) and the results concerning 
the random variable £ p . However, the increments of the process lnX p being 
independent, the convergence of its restrictions (and hence of those of X p ) 
on finite intervals [A, B] C R (that is, convergence in the Skorohod space 
T>[A, B] of functions on [A, B] without discontinuities of the second kind) 
follows from Gihman and Skorohod [H Theorem 6.5.5], Lemma [5] and the 
following lemma. 

Lemma 5 For any e > we have 

lim lim sup P< \\nX p (yi) — \nX p (y 2 ) \ > e \ = 0. 

/WO P^O \ yi -y 2 \ <h I ' J 

Now, Theorem [T] follows from the following estimate on the tails of the 
process X p by standard argument. 
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Lemma 6 For any b G ] , 3/40 [ we have 

p( sup X p (y) > e- bA \ ^ 2e~ bA 
l\y\>A ) 

for all sufficiently small p and all A > 0. 

All the above lemmas will be proved in the next section, but before let 
us discuss the second possible asymptotics p — > +oo. One can show that 
in this case, the process Z p converges weakly in the space T> (— oo, +oo) to 
the process Zqo^u) = e~ u t{ u>v j, u 6 R, where r\ is a negative exponential 
random variable with P{^ < t} = e*, t ^ 0. So, the random variables ( p and 
£ p converge weakly to the random variables 

L u Zooiu) du 

Coo = r ry / \ i =V + 1 and ^oo = argsup ZgoiU) = T) 
J R IvoKU) du ueR 

respectively. One can equally show that, moreover, for any k > we have 

ECj-EC^ and Eg - 

and in particular, denoting = E(^, = E^ and E^ = B^/M^, 
we finally have 5^5^= E(r/ + l) 2 = 1, M p -> = Er/ 2 = 2 and 
E p — > -Eoo= 1/2. 

Let us note that these convergences are natural, since the process Z^ 
can be considered as a particular case of the process Z p with p = +oo if one 
admits the convention +oo • = 0. 

Note also that the process Z^ (up to a linear time change) is the limiting 
likelihood ratio process of Model 1 (Model 2) in the situation when a ■ b = 
(S_ ■ S + =0). In this case, the variables Coo = V + 1 and £oa = V ( U P to a 
multiplicative constant) are the limiting distributions of the Bayesian esti- 
mators and of the maximum likelihood estimator respectively. In particular, 
Boo = 1 and M ra = 2 (up to the square of the above multiplicative constant) 
are the limiting variances of these estimators, and the Bayesian estimators 
being asymptotically efficient, E^ = 1/2 is the asymptotic efficiency of the 
maximum likelihood estimator. 

3 Proofs of the lemmas 

First we prove Lemma [2j Note that the restrictions of the process In X p (as 
well as those of the process In Zq) on R + and on R_ are mutually independent 
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processes with stationary and independent increments. So, to obtain the 
convergence of all the finite-dimensional distributions, it is sufficient to show 
the convergence of one-dimensional distributions only, that is, 

\nX p (y) =► \nZ (y) = W(y) - M = A^(-M , \y\ 

for all y £ R. Here and in the sequel "=^" denotes the weak convergence 
of the random variables, and Af(m, V) denotes a "generic" random variable 
distributed according to the normal law with mean m and variance V. 

Let y > 0. Then, noting that n + ^— J is a Poisson random variable of 

parameter A = — — ^ — > oo, we have 

P(e'-l) 

i n x pW = ,n + f^-^p 1 /=f=:5i^ + a » 



p/ p r Vp(e"-l) e"-l p 

- , r n +(f)~ A e'-l-p , Kr ( V 



since 



C-l >/A p(e^-l) 
p e^-l-p p 2 /2 + o(p 2 ) 1 



6 p_1 p + (p) ' p{eP-l) p{p + o(p)) 2 

and 

w -/V(o,i). 



Similarly, for y < we have 



lnX p (y) = -pIT_(^)-^ = p 



y\ y / -y -0 2/ 



p y p 'V P(l -e- p ) VA' l-e~" p 



and so, Lemma [5] is proved. 

Now we turn to the proof of Lemma H] (we will prove Lemma [3] just after). 
For y > we can write 

«^,-.- P (K ; )-i)-p(-i).-p(fn t (! 
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Note that n + ( — ) is a Poisson random variable of parameter A 



p> r />(<" i: 

with moment generating function M{t) = exp(A (e* — 1)). So, we get 



M v» to) = ^(_^^(_JL_ ( ^_ 1) ) 

= ex P -777 + 773775 , ^ = ex P 



3 P/2 _ ! 



2p p(e^/ 2 + l)7 2p(eP/ 2 + l 

e p/4_ e -p/4 \ / tanh(p/4)\ 

GXP 1 ^ 2p(eP/* + e-»/*) ) = GXP [~ y —Ip— ) 
For y < we obtain similarly 



E^(„) = E«p(-£n.(^)-X) 



i V , V \ ( 1 - ?-<" 2 



/ tanh(p/4)\ 

v 2p~; ' 



= exp 

Thus, for all y G R we have 



EA-/%) = exp(-|y|^|^) 



(5) 



and since 

tanh(p/4) _ p/4 + o(p) 1 
2p ~ 2p" * 8 

as p — > 0, for any c G ] , 1/8 [ we have EJfp ($/) ^ exp(— c\y\) for all 
sufficiently small p and all £/ G R. Lemma H] is proved. 



Further we verify Lemma [3j We first consider the case yi,y2 G R+ (say 
Vi ^ 2/2)- Using (j5J) and taking into account the stationarity and the inde- 
pendence of the increments of the process lnX p on R + , we can write 

E \X l J\ yx ) - Xy 2 ( y2 )\ 2 = EX p ( yi ) + EX p (y 2 ) - 2BX 1 / 2 ( yi )X 1 p % 2 ) 



2-2EX p (|/ 2 )E: 



Xl /2 (y2) 
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2-2EXy 2 ( yi -y 2 ) 
2 - 2exp ( - \yx - y 2 \ 



tanh(p/4)\ 



) 



P 4 

The case y\,y 2 G R- can be treated similarly. 
Finally, if y±y 2 ^ (say yi ^ ^ yi), we have 

e ix^feo - x;' 2 (y 2 )| 2 = 2-2 Ejfffe) ex^m 

= 2-2exp(H^-^f A 
1 . 

^ - \V\-Vt\ , 



) 



4 

and so, Lemma [3] is proved. 

Now let us check Lemma El First let yi, y 2 G R + (say yi ^ y 2 ) such that 
A = |yi - y 2 \ < h. Then 

p{|lnX,( yi )-lnX,(jfc)| >e} ^^EllnX^O-lnX,^)! 2 

lE|lnX p (A)| 2 



£ 2 



A\ A 2 



P ' P 



= l(p 2 (A + A 2 ) + ^-2AA 

= i (/5(P) A + 7 (P) A 2 ) 
< i (/3(p) + t(p) ^ 2 ) , 



where A = p ^_d is the parameter of the Poisson random variable H+ (— ) , 

/3(P) 



p(eP-l) 

P P 



(eP _i) p + (p) 
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and 



1 1 

7(P) = T~ n Tv?. + ~3. 



(eP-1) 2 p 2 p(eP-l) \p eP-l 
e p -l-p V ( p 2 /2 + o(p 2 ) \ 2 1 
p(eP-l)J \p(p + o(p))J "4 

as p — > 0. So, we have 

lim sup pjllnX^) -lnX p (y 2 )| > 5} ^ lim {P(p) h + 7 (p) /i 2 ) 
^° |yi-«2|<h 1 J p->-0 £ 2 



1 A /l 2 



-2 



and hence 



lim lim sup p(|lnX p (j/i) - lnX p (y 2 )| > e} = 0, 

/l-<-0 p^O |,„_„„ 1^1, I J 



1-2/2 |<h 

where the supremum is taken only over 2/1,2/2 G IR+- 

For 2/1,2/2 £ K_ such that A = I2/1 — 2/2I < ^ one can obtain similarly 

p{\\nX p ( yi )-\nX p (y 2 )\ > e] ^ ^E\\nX p ( yi ) - \nX p (y 2 )\ 2 

= j 2 {P'(p) A + 7'(P) A 2 ) 



< \ 2 {P\P) h + i{p) h 2 ) , 



where 

P P 



(1-e-P) p + o(p) 

and 



7'(P) 



e-P-l + p V _ / p 2 /2 + o(p 2 ) y 1 



as p — > 0, which will yield the same conclusion as above, but with the supre- 
mum taken over 2/1,2/2 £ 

Finally, for 2/12/2 ^ (say 2/2^0^2/1) sucn that 1 2/1 — 2/2 1 < h, using the 
elementary inequality (a — 6) 2 ^ 2 (a 2 + 6 2 ) we get 

p{|lnX p (2/i)-lnX p (2/ 2 )| >e} < ^ E|lnX p ( yi ) - lnX p (y 2 )\ 2 
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2 

< — 

" e 2 



(E\lnX p ( yi )\ 2 + E\\nX p (y 2 ) 



1 {Mih.^{p)&+P(p)\to\H(p)\v*?) 



e 2 
2 



which again will yield the desired conclusion. Lemma [5] is proved. 

It remains to verify Lemma El Clearly, 
p( sup X p (y) > e~ bA ) ^ p(su P X p (y) > e- bA ) + p( sup X p {y) > e~ bA 

l\y\>A ) ly>A ) ly<-A 

In order to estimate the first term, we need two auxiliary results. 

Lemma 7 For any c G ] , 3/32 [ we have 

EX p 1 / 4 (y)^exp(-c|y|) 
for all sufficiently small p and ally G R. 

Lemma 8 For all p > the random variable 

r] p = sup (lT A (t) - t), 

teR+ 

where U\ is a Poisson process on R + with intensity X = p/(e p — 1) G ]0, 1[, 
verifies 

Eexp (^ p ) ^2. 

The first result can be easily obtained following the proof of Lemma HI 
so we prove the second one only. For this, let us remind that according to 
Shorack and Wellner [191 Proposition 1 on page 392] (see also Pyke [17]) . 
the distribution function F p {x) = P{r] p < x} of rj p is given by 



1 - F p (x) = P{rj p > x} = (1 - A) e x * (A e~ x ) 



n>x 

for x > 0, and is zero for x ^ 0. Hence, for x > we have 

l-f»<(l-A)^E <pf J Xe-r 

~* v 27m n n e n 

n>x v 

14 



27T ~ y/n V 71/ 



2tt ^ V™ 



^l^ e (A-l>( Ae l-A)^( Ael A )' 



27r ~ \Jn- x 



1_ A ™ (Ae 1 -*)* 1-A f (Ae 1 "*) 1 
A" V ^ 1 J- < -=■ A" / ^ dt 



k>0 



i-x xx r(l/2) l-A 



2yr ^/-lnfAe 1 -^ J _21n ( Ael " A ) 
V2sinh(p/2) / 



- 1 



o P /2 _ p-p/2 



where we used Stirling inequality and the inequality 1 — A ^ y — 2 In (A e 1 A ) , 

which is easily reduced to the elementary inequality ln(l — fi) ^ — /i — /x 2 /2 
by putting \i — 1 — A. So, we can finish the proof of Lemma [S] by writing 

Eexp(^ p ) = / e^dF p (x) 

J [R 

= e^ /4 (F p (x) - 1)1 + °° - £ / e px/ \F p (x) - l) cfe 
L J -oo 4 Jj, 



/ e px/4 + % f e px/4 (l - F p {x)) dx 



< 1 + ~ / e^ px/4 cfe = 2. 

Now, let us get back to the proof of Lemma El Using Lemma El and taking 
into account the stationarity and the independence of the increments of the 
process lnX p on R + , we obtain 

p(supX p (y) > e- bA \ ^ e bA/4 EsuvXy\y) 

{y>A J y>A 

xl%) 



e^EX^(A)Esup 

y>A X p ' (A) 

e bA/A EX 1 / 4 (A) EsupXy\z) 

z>0 
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= e hA /' EXy*(A) Esup [e W ^U + (z/p) - J)) 

= e M / 4 EXy 4 (A) Eexp (sup (^(n - t))) 
= e M/4 EX 1 / 4 ^) Eexp r/ p ) < 2e bA/4 EX^A). 

Hence, taking b G ]0,3/40[, we have 56/4 G ]0,3/32[ and, using 
Lemma d we finally get 



P jsupXpfo) > e- 6A | ^ 2e 6A/4 exp(- = 2 



for all sufficiently small p and all A > 0, and so the first term is estimated. 

The second term can be estimated in the same way, if we show that for 
all p > the random variable 

i = sup (-II A ,(t) + t) = - inf (Uy(t) - t), 

teR+ teR + 

where Uy is a Poisson process on R + with intensity A' = p/(l — e~ p ) G ]0, 1[, 
verifies 

Eexp (|V P ) <2. 

For this, let us remind that according to Pyke [T7] (see also Cramer [2]), 
?7 p is an exponential random variable with parameter r, where r is the unique 
positive solution of the equation 

A'(e~ r - 1) + r = 0. 

In our case, this equation becomes 

P 



1 - e~P 



(e- r -1 +r = 0, 



and r = p is clearly its solution. Hence r]' p is an exponential random variable 
with parameter p, which yields 



Eex p(ii) 

and so, Lemma [6] is proved. 



4 

3 <2 ' 
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4 Numerical simulations 



In this section we present some numerical simulations of the quantities B pi M p 
and E p for p e ]0, oo[. Besides giving approximate values of these quantities, 
the simulation results illustrate both the asymptotics 

B M 
Bp ~ — , M p ~ — and E p -> £ as p -> 0, 
p P 

with 5 = 16 C(3) « 19.2329, M = 26 and £ = 8C(3)/13 « 0.7397, and 
Bp -> Boo, M p -> Moo and E p ~^ E^ as p -> oo, 

with .Boo = 1, Moo = 2 and = 0.5. 

First, we simulate the events xi,x 2 , ... of the Poisson process 11+ (with 
the intensity l/(e p — 1)) , and the events x[, x' 2 , . . . of the Poisson process n_ 
(with the intensity 1/(1 — e~ p )). 

Then we calculate 

xZp[x) dx 



Zp(x) dx 



' R 

oo oo oo oo 

p—pt+x'. 



1 i=l i=l j=] 

oo oo 

p—pi+x'- 



and 



i=i i=i 



Xfc, if pk - x k > p - p£ + a^, 



where 



£ p = argsup Z p (a;) = , 

xgr I — x e , otherwise, 



k = argmax (pi — Xj) and £ = argmax (p — pi + x^), 
i>i i>i ' 



so that 

Xfc = argsup Z p (x) and — = argsup Z p (x). 

x£R + xeR- 

Finally, repeating these simulations 10 7 times (for each value of p), we 
approximate B p = and M p = by the empirical second moments, 
and Ep = B p /M p by their ratio. 
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The results of the numerical simulations are presented in Figures [T] and [21 
The p — ► asymptotics of B p and M p can be observed in Figure [TJ where 
besides these functions we also plotted the functions p 2 B p and p 2 M p , making 
apparent the constants Bq rs 19.2329 and Mq = 26. 




Figure 1: B p and M p (p — > asymptotics) 



In Figure [2] we use a different scale on the vertical axis to better illustrate 
the p — > oo asymptotics of _B p and M p , as well as both the asymptotics of E p . 
Note that the function E p appear to be decreasing, so we can conjecture that 
bigger is p, smaller is the efficiency of the maximum likelihood estimator, 
and so, this efficiency is always between Eqo = 0.5 and Eq ~ 0.7397. 




Figure 2: B p and M p (p — ► oo asymptotics) 
E p (both asymptotics) 
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